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Abstract 



Q^, Effects of CP-violating phases in active-sterile solar neutrino oscillations are 

o,| discussed in a general scheme of 3+A^s mixing, without any constraint on the 

D ' mixing between the three active and the Ng sterile neutrinos, assuming only a 

"^ - realistic hierarchy of neutrino mass-squared differences. A generalized Parke for- 

mula describing the neutrino oscillation probabilities inside the Sun is calculated. 
The validity of the analytical calculation and the probability variation due to the 
f-^ . unknown CP-violating phases are illustrated with a numerical calculation of the 

^^ I evolution equation in the case of 3+1 neutrino mixing. 
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1 Introduction 

Despite the success of standard three-neutrino oscillations [1] in explaining the results of 
solar (SOL), atmospheric (ATM), reactor and accelerator neutrino experiments with two 
distinct mass-squared differences (i.e., ArrigoL and Am\rj,-^) and three non-zero mixing 
angles [J( i.e., 612, 623 and 6*13), some anomalies in short baseline (SBL) neutrino oscilla- 
tion experiments (e.g., LSND |2], MiniBooNE [3], the Reactor anomaly [1] and Gallium 
anomaly j5]) indicate the existence of oscillations with much shorter baselines. This would 
imply the existence of extra mass-squared differences (i.e., Arngg^) with the hierarchy 

and therefore the mixing of three active neutrinos with extra sterile neutrino states [6]- 
|8]. Furthermore, the analysis [HIITO] of cosmic microwave background and large scale 
structure data may hint at the existence of additional radiation in the Universe, with 
sterile neutrinos being one of the plausible candidates. One extra sterile neutrino is also 
allowed by recent analyses of Big Bang Nucleosynthesis [TT]. Therefore, we should be 
open-minded on the existence of sterile neutrinos and it might be instructive to study the 
effects of the light sterile neutrino hypothesis in solar neutrino p^lTSj and atmospheric 
neutrino [H] oscillation experiments or in non-oscillation processes including beta decay 
[T^ and neutrinoless double-beta decay [TB] . 

Solar neutrinos produced in the core of the Sun can undergo matter-enhanced Mikheev- 
Smirnov-Wolfenstein (i.e., MSW [TTtlTB]) oscillations when they propagate inside the Sun. 
When a generic scheme with one or more sterile neutrinos is considered, some sub-leading 
effects including the effects of light sterile neutrinos [19], non-unitarity of the lepton mix- 
ing matrix (NU) [20] and non-standard interactions (NSI) [21] may contribute to the neu- 
trino oscillation probabilities and modify the standard MSW picture in the three-neutrino 
mixing scheme. Therefore, it is possible to constrain or measure these high-order effects 
with future precision solar neutrino experiments. 

In an earlier work [12] by two of the present authors, we proposed a general method 
to calculate the matter effects of solar neutrino oscillations with an arbitrary number of 
sterile neutrinos, without any constraint on the magnitudes of the active-sterile mixing. 
However, the oscillation probabilities derived in Ref. [12] are only valid for a real neutrino 
mixing matrix because of an incorrect treatment of the CP- violating phases. As will be 
explained in Section |3l in Ref. [12] the effects of the CP-violating phases have been 
removed with an inappropriate phase-transformation. 

In this work, we study the evolution of the neutrino flavor amplitudes inside the Sun by 
taking into account the roles of the CP-violating phases of the neutrino mixing matrix. 
We calculate the generalized Parke formula describing the electron neutrino survival 
probability and the electron-to-sterile neutrino transition probability. Furthermore, we 
validate our analytical calculations through a numerical solution of the neutrino evolution 
equation in the case of four-neutrino mixing and we illustrate numerically the effects on 
the survival and transition probabilities of the three CP-violating phases. 

This paper is organized as follows. We review the general framework of neutrino flavor 
evolution in Sec. [2] and we present the analytical expressions for the neutrino oscillation 



^The mass-squared differences and the three mixing angles are defined according to the standard 
parametrization in the latest PDG publication [T]. 



probabilities in Sec. [31 In Sec. H] a numerical validation of the analytical results and an 
illustration of the effects of the CP-violating phases is presented in the simplest case of 
four-neutrino mixing. Finally, we conclude in Sec. O At the end of the paper there are 
three appendices on the analytical derivation of the non-adiabatic crossing probability, 
on the explicit parametrization of the mixing matrix and on the density matrix method. 

2 General Framework 

Following the notation of Ref. [12], in this Section we shall give a brief review on the 
general framework of neutrino flavor evolution with three active and Ns sterile neutrinos. 
The flavor eigenstates of active neutrinos and sterile neutrinos can be written as 

N 

K) = J2^:^W^), (2) 

i=l 

where A^ = 3 + Ns, a runs over e, /i and r for three active neutrinos and Si, . . . , s^^ for 
N^ sterile neutrinos, u- is one of the N mass eigenstates with mass rrii, and U^^ stands 
for an element of the (3 + N^) x (3 + N^) neutrino mixing matrix. In the general case, 
solar neutrinos are described by the state 

Wi^)) = Yl ^a(^)ka) ) (3) 

o=e,/^,T,si,...,sjVs 

where x is the distance from the production point during the propagation with the initial 
condition ipa{^) = ^ae and the normalization ^crl^"(^)l^ ~ -'-• "^^^ Mikheev-Smirnov- 
Wolfenstein (MSW) equation [T71[IB] describing the evolution of the flavor transition 
amplitudes ifja^x) is given by (see Ref. [22] ) 

i A^ = Uy^ = {UMW + V) ^ , (4) 



dx 



with 



^V 2E ' 2E ' 2E ' ' 2E y ' ^ ^ 

V = diag(V^cc + ^Nc, V^c, V^c, 0, . . . , 0) , (7) 

where E is the neutrino energy and ^fril- = m^ — m^ . The charge- current and neutral- 
current matter potentials are defined as 

ycc = V2GfN, - 7.63 X 10"^^ ^' eV , H^c = -^v^GpiV^ , (8) 

Aa cm~'^ 2 

where Gp is the Fermi constant, N^ is the electron number density, Nn is the neutron 
number density, and A^a is the Avogadro's number. By using the definition of electron 
fraction Y^, = Ne/{Ne + Nn), we have 

N, = ^NaY, , Nn = ^Nj,{l-Y,) , (9) 

g g 



in an electro-neutral medium, with p being the mass density. Thus, we have the following 
relation between the two matter potentials: 

V^c = -RncVcc , with i?Nc = — ^rr^ ■ (10) 

Next, we can introduce the vacuum mass basis as 



/V _ fj,y „/,v\^ 



vP^ =«,..., z^5;)^=f/tv^ , (11) 



which satisfies the following evolution equation 

d 
dx 



t—^^ = {M + U^VU) ^^ . (12) 



Then we can decouple the fiavor transitions generated by Arn^i from those generated by 
the larger mass-squared differences, according to the following hierarchjo 

^cc - iV^NcI - ^ « ^1^ for k>3 , (13) 

for the solar matter density. Therefore, the N-component evolution equation Eq. (IT^ 
can be truncated to 






(14) 



and 



V^r(x)=^V^r(0)exp(-z^^) , for A:>3 . (15) 



^ + ^E(l^-l' + |f^-n^" ' (16) 



By subtracting a diagonal term 

AE 

a 

which generates an irrelevant common phase, the evolution equation in Eq. flT^ can be 
written as 

^^^y = HM , (17) 

with^y = (V'y,^y)^ and 

^■^ ~ • r sin 2^6-*"^ (5-1/cos2e 
where the variables 6, V and ^ are defined as 



Y _ 1-6 + Vcos2^ V sin 2^e'^ 



^ The different case of active-sterile neutrino mixing with a much smaher mass-squared difference in 
solar neutrino oscillations, has been studied in Ref. ^^. 



with 



Vcos2^ = J2{\Uai\^-\Ua2\^)^c. = lvccX , (20) 

a 

Vsm2^e'^ = Y,U:iUa2Vo^ = ^VccY = ^Vcc\Y\e'^ , (21) 

a 

X=|f/elp-|f/e2|' + i?NC Yl {\Ual\' ' \U^2\') 

Ns 
= \Uel\^ - |f/e2r - RnC J^ (l^^^ll' " l^-^^l') , (22) 

= 2 (u:,U,, - i?Nc f; KiU,^,] . (23) 



«=1 



Therefore, we can obtain the full expressions for ip, C, and V as 



<^ = arg(F) , (24) 

tan2e = ^ , (25) 

V = WccVx^T]yf . (26) 



To be more explicit, we can rewrite the Hamiltonian in Eq. f lT8|) in a compact form: 

n^ = M2 + UIV2U2 , (27) 

with M2 = diag(— 5, 5), V2 = diag(V, — V^) and 

U2 = W2{^,^)^( ^°^/_., ^^"^f) , (28) 

where 1^2(^5 V^) is the complex rotation matrix which can be generated by a real rotation 
matrix R2{0 ^^^ ^ diagonal phase matrix D2{(p) (see Ref. [22]), 

W2{^,^) = Dl{^)R2{0D2{^) , (29) 

i?.(0-f'''S "'fl , (30) 

^^ \^— sm4 cost, J 

D2{v)=dia.g{l,e'n . (31) 

Using the definition of these effective parameters, we have obtained an evolution equation 
analogous to that in the two-neutrino mixing scheme. But one should keep in mind that 
the effective mixing angle ^ and phase (p are not simple mixing parameters given by 
a specific parametrization but medium- dependent parameters which are constant only 
when the electron fraction Y^ remains unchanged along the neutrino propagation path. 



To solve the evolution equation in Eq. ( 1T7|) . we can first diagonalize the Hamiltonian 
in Eq. fl27|) with a complex rotation 



m^^=W2{u,^)m^ , (32) 

where 

ysin2e .„„. 

— K cos 24 

(f is the complex phase defined in Eq. (1241) . and \l/^ is the amplitude vector in the effective 
mass basis in matter. Then the evolution equation becomes 

where the Hamiltonian can be decomposed into the adiabatic (ad) and non-adiabatic (na) 
parts, 

njM _ njM , njM 
'^2 — 'T-ad "T ''-na 

6mJ \{l^sm2uj + iuj)e-'^ ^sin^w J ' ^ '' 

with 



5m = V(5 - V cos 2^)2 + (y sin 2^)^ 

= V(y-(5cos202 + (5sin202 , (36) 

. _dip . _duj 



(3^ 



Finally, we can arrive at a formal solution of Eq. ([3; 



|^2^(X,)|2 = K(0)pPi2 + li^fm' (1 - P12) 



for the averaged amplitudes of solar neutrino evolution, where Xd is the coordinate of the 
detector on the Earth, P12 is the level-crossing probability between two effective mass 
eigenstates ip^, ip^ during their propagation inside the Sun. By definition, P12 is zero 
when the off-diagonal terms in Eq. (135|) are vanishing, which is defined as the adiabatic 
approximation. 

Before finishing this Section, we want to point out that, as an improvement to the 
results in Ref . [12] , in this paper we consider the general case in which the effective phase 
if cannot be absorbed by a simple rephasing transformation when the electron fraction 
Ye is not constant along the neutrino path. In this case, the CP-violating phases in the 
mixing matrix may influence the neutrino flavor evolution inside the Sun through the 
effective phase ip, which is the complex argument of Y Eq. fl2^ . Moreover, since also the 
module of Y depends on the CP-violating phases in f/^i and Ua2, these phases can affect 
also on the behavior of ^, V and uj. In the next Section we derive the analytic expression 
for the average oscillation probabilities of solar neutrinos taking into account the effects 
of the CP-violating phases in the mixing matrix. 
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3 Generalized Parke Formula 

In this Section, we derive the neutrino oscillation probabilities based on the framework 
presented in the previous Section. Due to the energy resolution of the detector and the 
uncertainty of the production region, the interference terms between the massive neutrinos 
are not measurable [23]. As a result, we obtain the averaged oscillation probabilities 



p: 





N 

k=l 


2 


N 




k=i 

2 N 


VflV^.^W ' 


fc=l 


fc=3 



(39) 

where the matter effects in the detector are neglected. According to Eq. flTSj) . ^^ with 
(A; > 3) is decoupled from other amplitudes inside the Sun and evolves independently. All 
we need is to solve the evolution equation in the truncated 1-2 sector, given in Eq. flT7|l 
in the vacuum basis or in Eq. ( 134|) in the effective mass basis in matter. From the initial 
conditions for the flavor amplitudes ipi3{0) = 5pei one obtain that in the vacuum basis 



ek' 

Writing Uj3i and Up2 as 

TT ■ a i<t>R9 with COS Xl3 = \Ul3l\ + f//32 , (40) 

(7/32 =sm0i3 cosx/3e^' v.^ 1 m 1 1 ^ 1 ^ v / 

the initial conditions in the effective mass basis in matter are 



^M(O) = cosXe(cosu;°cos6'ee"^'^'=i - sina;°sin6'ee-*('^<=2-v'') 
iPf{0) = cosXe(sina;°cos6'ee-*('^°+'^'=i) + cos u;° sin 6^66 



Upe2 ^ 



(41) 



where 0;° and y?" are the rotation parameters [see the definition in Eq. (132|) ] between the 
vacuum mass basis and effective mass basis at the production point. 

Using the relations in Eq. ( 138|) . we obtain the averaged solar neutrino oscillation 
probabilities 

N 

75^ 2 2 75(2^') I Y^lrr i2itt |2 / An\ 

Py.^v, = COS XeCOS X^Py^-^u. + } \Uek\\Upk\ , (42) 



fc=3 



with 



P^uliup = ^ + (^ - ^12) COS 2^/3 [cos 26 e cos 2a;° - cos $° sin 26 e sin 2a;°] , (43) 

where 

$° = 0el - 0e2 + / • (44) 

One can check that the survival probability reduces to the well-known Parke formula 
[25] in the limit of two-neutrino mixing, in which cos^ Xe = cos^ X/3 = cos $° = 1 and 



[cos 29e COS 2cij° — COS $° sin 29e sin 2a;°] — ?■ cos 2^°, where 6'° = O^+u^ is the effective mixing 
angle between the ffavor basis and the effective mass basis in matter at the production 
point. 

The oscillation probabilities reduce to the ones in Ref. [12] in the case of CP invariance, 
in which cos<l>° = 1. Note, however, that in the case of CP violation the expression fH2]) 
for the oscillation probabilities do not reduce to the corresponding Eq. (61) of Ref. [12] 
even if the electron fraction Ye is constant along the neutrino path. In this case, as 
explained in Ref. [12] one can eliminate the phase ip by an appropriate rephasing of the 
mixing matrix. In fact, since the MSW evolution equation (|4]) is invariant under the 
phase transformation 

U^k^e'^P-Ue^ue'^P^ (45) 

which transforms ip — ?■ ip + ip2 — <^i., one can eliminate if by choosing ^pi — ip2 = 'P- 

However, in this case there is no remaining phase freedom to eliminate both 0ei and 0e2, 

contrary to what has been incorrectly stated in Ref. [12]. Indeed, since 0ei ~ 0e2 + V' is 

g 

invariant under the phase transformation (1451) (as well as X and |y|, and hence P^^^y )■, 

it is clear that it cannot be eliminated. Hence, the oscillation probabilities in Eq. fH2|) 

represents the improvement of the corresponding oscillation probabilities in Eq. (61) of 

Ref. [12] which takes into account in a proper way the effect of the CP-violating phases 

in the mixing matrix |2|. 

The CP-violating phases can contribute to the oscillation probabilities in several dif- 
ferent aspects. It is obvious to identify cos$° as a term which gives a direct impact, but 
in practice cos ^^ is very close to unity for solar neutrinos, as shown in the next Section 
(see, for example. Fig. [7]). However, there can be a significant phase dependence of the os- 
cillation probabilities coming from a;°, which depends on the CP-violating phases through 
X and \Y\ and through the modules of the mixing matrix elements in two different rows, 
which depend on the cosines of the CP-violating phases. 

When counting the relevant number of CP-violating phases in solar neutrino oscilla- 
tions, one should take into account that v^ and v-r are indistinguishable and all sterile 
neutrinos are indistinguishable. Therefore, if the mixing matrix is written as a prod- 
uct of complex rotations, the rotation in the t'^-^'r sector and all the complex rotations 
among sterile neutrinos do not have any effect on the observables in solar neutrinos. 
These complex rotations can be eliminated from the evolution equation (TT2l) by choosing 
a parameterization of the mixing matrix in which they occupy the left-most positions. 
Then, the remaining part of the mixing matrix can be written in terms of two complex 
rotations among active neutrinos and three complex rotations between the three active 
neutrinos and each sterile neutrino, with a total of 2 -|- 3A^s complex rotations. Of the 
corresponding phases, there are 2 + Ng Majorana phases which can be factorized in a 
diagonal matrix on the right of the mixing matrix (see Section 6.7.3 of Ref. [22]) and 
have no effect on oscillations. Therefore, solar neutrino oscillations depend on 2Ns Dirac 
CP-violating phases. 

In conclusion of this Section, let us discuss the problem of calculation of the crossing 
probability Pi2- In Appendix |A] we describe two different approximations which allow us 
to derive analytical expressions for P12. The first method [12], described in Appendix lA.lt 



•^Some effects of the CP-violating phases in solar neutrino active-sterile oscillations have been discussed 
earlier in Refs. [5111^. 
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Figure 1: The magnitudes of four different terms in Eq. (135|) : Ti = 5yi, T2 = 10 sin ^ 



a; 



T3 = I i0 sin 2a; I and T4 



\0J\ 



The mixing parameters are set to Ml and PI in Ap- 



pendix |B] and the matter density distribution is taken from the BSB2005(OP) Standard 
Solar Model |29]. 



can be applied when the electron fraction Y^ is approximately constant along the neutrino 
propagation path. In this case, we can employ the similarity of the neutrino flavor 
evolution in Eq. fl34p to that of two-neutrino mixing and obtain the crossing probability 
with the help of non-perturbative calculations [2B]. The second method, described in 



Appendix \K.2\ can be applied when the non-adiabatic contribution in the Hamiltonian 
fl35|) is much smaller than the adiabatic one. In this case, one can use the general 
perturbation theory [271128] and calculate the effective crossing probability to include 
the non-adiabatic contribution along the whole path of neutrino propagation inside the 
Sun. 

Both of the above approximations are quite good for the solar neutrino evolution inside 
the Sun. As shown in Fig. 1 of Ref. [12], the electron fraction Ye is almost constant in 
the radiative (0.25 ^ r < 0.7) and convective (r > 0.7) zones, with the only exception of 
the core region with r < 0.25. However, the effects of the variation of Yg in the core are 
negligible if the flavor transitions occur mainly in a resonance located in the radiative or 
convective zone. The numerical analysis in Ref. [12] validated this approximation. For 
the approximation of perturbative expansion, we show the magnitudes of different terms 
of Eq. fl35l) in Fig. [H considering, as an example, the values of the mixing parameters Ml 
and PI in Appendix [B] and using the matter density distribution in the BSB2005(OP) 
Standard Solar Model [22] . One can see that the order of magnitude of the non-adiabatic 
terms can only reach at most about 0.1% of the adiabatic term, verifying the accuracy of 
the perturbative approximation. In practice, the numerical calculations of P12 with the 
two methods are consistent and both show that P12 is negligibly small. Therefore, in the 
numerical analysis discussed in the following Section we neglect the crossing probability 
P12. 
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Figure 2: Energy spectra of the analytical and numerical evaluations of the solar neutrino 
electron survival (upper panel) and electron-to-sterile transition (lower panel) probabili- 
ties. All the oscillation parameters are set to Ml and PI in Appendix |Bl 

4 Numerical Discussion 



In this Section, we illustrate the validity of the neutrino oscillation probabilities in Eq. 
and the effects of the CP- violating phases by using a numerical calculation of the neutrino 
evolution equation in the scheme of four-neutrino mixing. The explicit parametrization 
of the 4x4 neutrino mixing matrix and the values of the oscillation parameters used in 
the discussion are presented in Appendix [Bl 

As explained in Section [31 in the case A^^ = 1 the observable effects of solar neutrino 
oscillations depend on two Dirac CP- violating phases. In fact, we can write the mixing 
matrix as U = W'^^U', where W"^^ = W^(^23,''723) is defined in Eq. (165|) and U' is a proper 
product of the other rotations, which contains two Dirac CP- violating phases. Since W"^^ 
drops out of the evolution equation ( !T2|) . the oscillation probabilities are independent 
from ?723 (as well as from ^23)- They depend only on the two Dirac CP- violating phases 
in U'. However, in the following we will discuss the possibility to reveal the effects of 
the phases in a future scenario in which the absolute values \Ua4:\ of the elements of the 
mixing matrix with a = e, fi,T have been determined by precision short-baseline neutrino 
oscillation experiments. Hence, we adopt the parametrization in Appendix [B] in which 
\Ue4\, \Ufj,4\ and \Uri\ are independent of the phases and determine the mixing angles 
6'i4, 6*24 and 634 (there is no way to get such result with U = W'^^U'). Hence, although 
in the following we consider the three CP-violating phases in the parametrization in 
Appendix [HI one should keep in mind that the oscillation probabilities depend only on 
two phases, which are complicated functions of the three CP-violating phases and of the 
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mixing angles in the parametrization in Appendix |Bl 

We employ the data of the matter density distribution in the BSB2005(OP) Standard 
Solar Model [29] and we consider, for simplicity, neutrinos produced at the solar center. To 
obtain the numerical evolution of solar neutrinos inside the Sun, we use the fourth-order 
Runge-Kutta method described in Numerical Recipes [HI]. Since the unitarity condition 
is not automatically guaranteed in a straightforward application of the evolution equation 
in Eq. (IT7I) and can be violated by the errors of the numerical computation, especially 
for the evolution in the crucial resonance region where the amplitudes oscillate rapidly, 
we employ the equivalent density matrix formalism (e.g., see Chapter 9 in Ref. [22]) in 
which the unitarity condition is fulfilled by definition. We refer to Appendix O for a brief 
introduction on the basics of the density matrix method. 

In Fig. [2] we compare the analytical forms (ANA) of the electron neutrino sur- 

g 

vival probability Pge = Py^^^^ and the electron-to-sterile neutrino transition probability 

g 

Pes = Py^^y^ given by Eq. fH2l) with the corresponding numerical evaluations (NUM) 
of the neutrino flavor transitions. The upper and lower panels represent the electron 
survival and electron-to-sterile transition probabilities, respectively. We illustrate the 
comparisons with solid lines and cross points for the analytical and numerical oscillation 
probabilities, which show a perfect agreement between two different calculations of the 
evolution equation. Numerically, the accuracy of the analytical calculation is better than 
10~^ and no systematic deviation appears. 

Next, we want to illustrate the effects of the CP- violating phases in solar neutrino 
active-sterile oscillations. We can observe from Eq. fj^2|) that the oscillation probabilities 
are only sensitive to the absolute values of the elements of the neutrino mixing matrix, 
except for the explicit contribution of the phases in cos^*^. However, as discussed in the 
last Section, the CP-violating phases determine also the contributions of the modules 
of the elements of the mixing matrix, since two distinct rows (i.e., the electron and 
sterile rows) are involved in the oscillation probabilities. In the specific parametrization 
of the neutrino mixing matrix presented in Appendix [HI the modules of the matrix 
elements in the electron row are independent of the CP-violating phases, but those in 
other rows are phase- dependent. Therefore, the effects of the CP-violating phases in the 
electron neutrino survival probability arise only in the effective two-neutrino oscillation 
probability in Eq. 0431) via the effective mixing parameters u^ and $°. On the other hand, 
the CP-violating phases manifest themselves in the electron-to-sterile neutrino transition 
probability by the phase dependence in u^ , $°, 6*^, Xs- Note that the constant term in 
Eq. (H2l) (i.e., ^^=3 It^efePl^sfeP) depends on the variations of the CP-violating phases 
via Xs in our specific parametrization. 

To show the variation of the oscillation probabilities for different values of the CP- 
violating phases, we can measure the possible size of the probability variation as the 
difference between the maximal (MAX) and minimal (MIN) values of the probabilities 
in the full parameter space of the three CP-violating phases. Therefore, we define the 
following asymmetries of the oscillation probabilities: 

MAx[g[-^MiN[g] 

^^^) = ' ^ MAX[g] + MIN[Q] ' ^^^^ 

where Q could be either the survival or transition probabilities. Each asymmetry illustrate 
the possible variation of the corresponding probability depending on the unknown values 
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Figure 3: Energy spectra of Pee(^) — Pee(O) (upper panel) and A{Pee) (lower panel) for 
the electron neutrino survival probability. The mass and mixing parameters are set to 
Ml in Appendix [B] and the three CP- violating phases are randomly scanned in the full 
parameter space. The short-dashed line in the upper panel corresponds to the set of 
phases PI in Appendix [Bl 




10 20 



E (MeV) 



Figure 4: Same as Fig. [3l but for the electron-to-sterile neutrino transition probability. 
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Figure 5: Energy spectrum of cos26g in Eq. (H7|) . The mass and mixing parameters are 
taken as Ml in Appendix [B] and the three CP- violating phases are randomly scanned in 
the full parameter space. 



of the CP- violating phases in a future scenario in which the absolute values \Ua4\ of the 
elements of the mixing matrix with a = e, fi,T have been determined by precision short- 
baseline neutrino oscillation experiments [19]. In the following numerical discussion we 
consider, as a realistic example, the values Ml in Eq. f l75]) of the mixing angles, which 
determine the absolute values of the relevant elements of the mixing matrix through 
Eqs. dHZD-lIZI]). 

In the lower panels of Figs. [3]and|4]we show the energy dependence of the asymmetries 
A{Pee) and A{Pes). In the upper panels we show the possible range of variation of the 
probabilities Pee and Pes for all possible values of the phases 7714, 7724, 7734, with respect 
to the case 7714 = 7/24 = 7/34 = 0. The shadowed regions with red and green colors are 
generated by scanning the full parameter space of three CP-violating phases. The two 
boundary curves stand for the maximal and minimal values of the differences (which 
correspond to the maximal and minimal values of the corresponding probability). These 
maximal and minimal values are used in calculating the corresponding asymmetry A{Pee) 
or A[Pf,s) in the lower panels of Figs. |2] and HI Note that the boundary curves may 
correspond to different values of the CP-violating phases for different energies. 

In the upper panels of Figs. [3] and H] we have also shown the curves corresponding 
to the values PI in Eq. (174|) of the CP-violating phases. We can observe that the vari- 
ation induced by these values of the three CP-violating phases is less than 1.2% for the 
electron survival probability and can be as large as 100% for the electron-to-sterile tran- 
sition probability. This is because the phase- independent contribution dominates in Pee, 
whereas both the phase-independent and phase-dependent contributions are comparable 
in Pes and both can induce significant variations in the transition probability. Notice that 
there is a kink and a sudden turn at about 2 MeV in the spectra of A (Pee) and A{Pes), 
respectively, which correspond to the similar behaviour of the MAX boundary curves in 
the upper panels. This property can be understood with the help of Fig. which shows 
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Figure 6: Scatter plots of the asymmetries AiP^e) and ^(Pe^) with respect to the con- 
tributions of the first two mass eigenstates (Pgg or P^^) and those of all the other mass 



eigenstates (cos^ Xe or cos^ Xs 



the energy spectrum of the quantity 



cos 29° = cos 2ee COS 2uj° - cos $° sin 2^^ sin 2uj° 



(47) 



in Eq. P3|) with randomly scanned CP- violating phases. Since cos 26° changes sign at 
about 2 MeV, the values of the CP- violating phases which maximize the probability have 
a sudden jump, which generates a sudden change of the slope of the curve of maximal 
probability. 

Let us now discuss separately the different contributions to the probability variation. 
In Fig. [6] we show the scatter plots of A{Pee) versus cos^ Xe and Pgg and of A{Pes) versus 
cos^ Xs and Pg^ obtained with a random generation of the three CP- violating phases in 
the entire parameter space. We considered a neutrino energy of 10 MeV and the mixing 
parameters Ml in Appendix [Bl Fig. [6] shows that the effects of the CP-violating phases 
show up only in the effective two- neutrino probability Pgg for the survival probability Pee, 
but emerge in both Pg^ and the suppression factor cos^ Xs induced by the other neutrino 
states for the transition probability Pes- Moreover, we can further study separately the 
phase dependence of Pgg and Pg^ due to the phase dependence of cos^6'e/s, cosa;° and 
cos$°, which is illustrated in the scatter plots in Fig. [71 In our parametrization of the 
mixing matrix, cos 26e is independent of the CP-violating phases, but cos 26s can reach 
almost all the possible values with the varying phase parameters. Therefore, the variation 
of the two-neutrino survival probability Pgg is dominated by the changing of cosa;°. 
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but the variation of the two-neutrino transition probabihty Pg^ conies from the phase 
dependence of cos 29s and cos w". On the other hand, from the rightmost panels of Fig. [TJ 
cos $° is very close to unity in both the survival and transition probabilities. This fact can 
be explained by considering Eq. ([21]), where the imaginary part of Y is suppressed by both 
the small active-sterile mixing and the small neutral current contribution (-Rnc — —0.2). 
Finally, the scatter plots in Fig. |8] show the correlations of the variations of cos u^ with 
those of X^ and \Y^\. One can see that the variation of cos uj^ is mainly determined by 
X^ rather than \Y^\, because of the small phase dependence of \Y\. 

In summary, we can conclude that the variation of the survival probability Pee due 
to the unknown CP-violating phases in the mixing matrix is determined mainly by the 
contribution of w'', whereas the transition probability Pes is sensitive to the CP-violating 
phases via Xsi (^s and uq, and the most significant contribution comes from Og. We 
have also shown that the direct phase dependence of the probabilities through cos $° in 
Eq. (H3|) is negligible because cos $° is very close to one in the full parameter space. 



5 Conclusion 

In this work, we have calculated the analytical solution of the flavor evolution of solar 
neutrinos in a general scheme of 3+A^s neutrino mixing, without any constraint on the 
mixing between the three active and the A^^ sterile neutrinos. We have improved the 
previous study in Ref. [12] by including the possible roles of the CP-violating phases 
in the mixing matrix and we have discussed the effects of these phases in active-sterile 
neutrino oscillations. We derived generalized Parke formulae which are suitable to be 
used in future precision measurements of solar neutrino oscillations. 

In Section Hj we have presented a numerical discussion with a realistic example of the 
possible phase contribution to the oscillation probabilities in the case of 3-1-1 neutrino 
mixing. We validated the analytical formulae with a careful numerical solution of the 
evolution equation inside the Sun. We illustrated the effects of the CP-violating phases 
through an appropriate asymmetry of the oscillation probabilities. We have shown that, in 
our example, the variations induced by the three unknown CP-violating phases can reach 
the level of 1% for the electron survival probability and may be as large as 100% for the 
electron-to-sterile transition probability. This scenario will be realized when the absolute 
values of the elements of the mixing matrix \UaA\ for a = e,fi,T will be measured in 
precision short-baseline neutrino oscillation experiments. In this case, it might be possible 
to observe the effects of the CP-violating phases in future solar neutrino experiments. 
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A Analytical Derivation of P] 



12 

In this Appendix, we present two methods for the approximate calculation of the crossing 
probability Pi2'- the constant Kg approximation and the approximation of perturbative 
expansion. 

A.l Constant Y^ Approximation 

Let us consider a case in which Yg is approximately constant inside the Sun. Then, -RnC) 
C, and if remain approximately unchanged during the neutrino propagation. Therefore, 
we can introduce the tilded vacuum mass basis defined by 

^y = Z}2M^y , (48) 

to accommodate the phase ip inside the amplitude vector {D2{ip) is defined in Eq. ( 131]) ). 
In the new basis, the evolution equation in Eq. (TT7|) becomes 

with§y = (^V^^V)T ^^^ 

^^-[ l^sin2e 6-Vcos2^) ' ^^^> 

Meanwhile, we can decompose the Hamiltonian into the vacuum and matter parts as 

n^ = M2 + UIV2U2 , (51) 

with U2 = i?2(0 [defined in Eq. (130|)]. 

As in the discussions in Ref. [12], we can further introduce the tilded effective inter- 
action basis \&2 ^^d the tilded effective mass basis \^^ in matter defined by 

¥, = R2i0^f, ^^ = R2{co)^f. (52) 

In this way, we obtain evolution equations for '$\ and \1/^ which have the same form as 
those without the CP- violating phases [12] if all the quantities with tildes are replaced by 
those without tildes. For instance, in the ^^ basis the evolution equation can be written 
as 

dx '' \ lUJ Sm J 

which is just the standard evolution equation of two- neutrino mixing in matter [22]- From 
the similarity we can define the adiabaticity parameter 

^ = ^ = '^ (54) 

' \uj\ ■l/5sin2^|dlniVe/dx| ' ^ ^ 
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and obtain the crossing probability Pu as 

exp (-f 7rF) - exp (^-f 7r ^ 



sirr^ 



1 - exp (-f 7R-^ 



9{Vo-Vn) , (55) 



where Vr = 5 cos 2,^ defines the resonance point and 7r is the adiabaticity parameter at 

the resonance with 

^ 2^sin^2e 

^"^ cos2e|dlniVe/dx|j^ ' ^ ' 

Finally, the 6 function is used to reduce P12 to zero when the potential at the production 
point is smaller than that at the resonance point. We use F = 1— tan^ ^ for an exponential 
density profile, which is a good approximation for the solar neutrinos 



A. 2 Perturbative Expansion 

As discussed in Section 3, the non-adiabatic terms are much smaller that the adiabatic 
term and we can treat Ti^g^ in Eq. fl5S]l as a perturbation term relative to "H^. Therefore, 
we can solve the S'-matrix defined in the effective mass basis, 

^f{xf) = S^{xf,x,)^f{x,) , (57) 

by using the standard perturbation theory (see Appendix B of Ref. [28]). After a straight- 
forward calculation, we arrive at the expression of 

^perti^f^^i) — ^a.d{^f^^i) ~ '^^a.di^f^^i) / '-'adl^'^i) '^na(^)'-'ad(^' ^V^^: 

Jxi 



' ' C* A 



^a.d{^ f^^ij ^^adi^f^^i) 



(1 + iA)e'^ -iCe'^ 
-iC*e-'^ (1 - iA)e-'^ 



(5^ 



where 



g ^-iA{xf,Xi) 

A(xf,Xi) = I 0sin^a;dx , (60) 



(59) 



C(x/,Xi) = /" ^(-0sin2a;-zc;;)e*(^-2^(^'^»)Mx , (61) 

Jxi 2 

A{xf,Xi) = / ^Mda; . (62) 



Xi 



Then, the effective crossing probability P12 is given by the probability of 1 ^ 2 non- 
adiabatic transitions 



Pi 



12 



{s^,A^f,o)} 



12 
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\C{xj,0)\' 



'M 



V5sin2a;-za;)e^(^-2^(^'°»dx 



(63) 



B Explicit Parametrization of U 

The 4x4 neutrino mixing matrix can be parametrized (see Ref. [22] for detailed discus- 
sion) as an extension of the standard parametrization [1] of three-neutrino mixing: 



34TI/24 d14 D23T;r/13 d12 



U = W^'^W'^R'^R'-'WR 



(64) 



where W°''^ = W{9ab, r]ab) and i?"* = W^^i^Oab, 0) are the complex and real unitary matrices 
in the (a, 6) plane, where W{9ab,f]ab) is defined by 



[W{6ab, r]ab)]rs = Srs + {cOS 6ab - l){6raSsa + Srh^sb) 

+ sm9abie-'''''''SraSsb - e'^^-^'drbSsa) 



(65) 



with 6ab and rjab being the mixing angles and Dirac CP phases in the specific plane. 

In this parametrization, the explicit expressions for the elements in the electron and 
sterile rows of U are given as follows 



Uel — C12C13C14 , 



Ue2 — ■S12C13C14 

Ue4 = Si4 , 



(66) 
(67) 



Usl — — S14C12C13C24C34 + (S12C23 + Si3e*''"S23Ci2)S24e*''^*C34 

+ (-S12S23 + Si3e*^^^Ci2C23)s34e*^^^ , 
Us2 = - S12S14C13C24C34 + (-C12C23 + 5125136*''" S23)s24e''''*C34 

+ (S23C12 + Si2Si3e*^"C23)s34e^''^^ , 
Us3 = - Sl3e"*''l^Si4C24C34 - (s34e*'''"C23 + 523^246^''"" 634)013 , 
Us4 =614624634 . 



Moreover, we have 



U, 



fj,4: — 6148246 



-«'?24 



f/r4 — ^146248346 



■IVM 



(68) 

(69) 
(70) 
(71) 



(72) 



In our numerical calculations, we consider the following values of the oscillation parame- 
ters: 



( Am 



Ml 



12 



7.54 X 10-5 eV 



^12 - 33.6° 
^23-39.r 
^13 ~ 9.0° 

^14 = ^24 = ^34 



(73) 



10°, 



and 



PI: r/i3 = 35°, r/24 = 75° , r/34 = 115° 



(74) 
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where the A?n^2 ^^^ three mixing angles (6*12, 6*13 and ^23), equivalent to the case of three- 
neutrino mixing, are taken from the latest global analysis [30], the active-sterile mixing 
angles are motivated by the anomalies of SBL data [2H5] and the phases are chosen 
non-trivially to reveal the effects of the CP phases. The assumed values of the active- 
sterile mixing angles 6*14, 624,, ^34 do not significantly affect the values of the oscillation 
parameters of active neutrinos extracted from the current data. 

C Density Matrix Method 

The density matrix formalism is equivalent to the framework of flavor amplitude evolution 
in Section 2. We consider that a neutrino state at the position x is described by the 
Hermitian density matrix operator 

p{x) = Y,Wa{x))W^Mx)\ , (75) 

a 

where Wa is the initial statistical weight of flavor a (i.e. the probability of the flavor a 
at X = 0). One can choose Wq, = 5ap for a neutrino state of one initial flavor /3. In the 
flavor basis we can define the specific density matrix as 

pI^ = {^v\p{^)\^d = ^We,'4)a^{x)i)l^{x) . (76) 

a 

The evolution equation of the density matrix p^ in the flavor basis, obtained from the 
evolution equation in Eq. (jlj), is 

■ P 1/F F Fl/F /^7^7^ 

^"d^ = ^ P ~ P ^ ' (77) 

with the initial condition pL(0) = Wr^6r)^. The density matrix in the vacuum mass basis 
p^ = U'^p^U follows the evolution equation 

,^=7/V^v_^v^v ^ (78) 

dx 

with 'Hy = U'^'H^U . For solar neutrino oscillations, using the approximation in Eq. f lT3|l . 
we can obtain the reduced evolution equation 

■^P2 1/V V Vi/V /nr\\ 

^-^ = ^2 P2 - P2 ^2 , (79) 

in the 2x2 subsystem of {ipY^ip^), where "H^ is defined in Eq. (TT7|) and p^ is given 
by {P2{x)}ij = ipY{x)ipJ (x) . Since "H^ is Hermitian and traceless, from the initial 
condition Wa = Sea we have 

Tr[py(x)] = iV'rr+iV'2''r = if/eir+it/e2r • m 

Therefore, the unitarity condition is fulfilled by definition. To be more explicit, we can 
rewrite these matrices in the terms of Pauli matrices with 

n^ = -^ayB , (81) 
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P2 = ^ 1 + -(TV ■ ^ . (82) 

where 

3 3 3 

?v = E^«^"'v' ^ = 5]i?:e«, S = Y,SjS^ , (83) 

a=l a=l a=l 

with (Ta (a = 1,2,3) being the Pauh matrices and (eyjeyjey) being three orthonormal 
vectors which form the vacuum mass basis. The components of the vectors B and S in 
the vacuum mass basis are 

B = (-2 Re{-Hy}i2, 2 lm{Hj}i2. {•H^}22 - {^^}ii ) (84) 

S={2 Re{py}i2, -2 Im{pJ}i2, {pI}ii - {pI}22) • (85) 

The evolution equation of the vector 5* is 

^ = 5^x5 , (86) 

ax 

with the initial condition 

5^(0) = (2Re(f/:i[/e2), -2lm{U:^U,2), Wei? - \Ue2?) . (87) 

According to Eq. f p9|) . the oscillation probabilities can be written as 

2 N 

Pl^u, = Y.\Upk?{pU^f)U + Y.\U^k?\Ue,? . (88) 

fc=l fc=3 

Using Eqs. fl86j) and fl88j) . we can employ the fourth-order Runge-Kutta method to perform 
the numerical evaluation of the neutrino flavor evolution. 
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